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ABSTRACT 


Adaptive  detection  and  copy  of  temporally  coherent  narrowband  waveforms 
with  multiple  sensors  is  discussed.  Unlike  related  problems  treated  in  the  radar 
literature,  the  sensor  response  to  the  coherent  signal  of  interest  is  assumed  unknown. 
Thus  the  techniques  presented  can  be  used  with  uncalibrated  sensor  arrays.  The 
signal  processing  takes  advantage  of  the  sensor  array  as  well  as  the  signal  coherence 
in  order  to  form  an  adaptive  spatial  beam  that  suppresses  cochannel  interference. 
More  general  signal  models  involving  partial  coherence  are  also  considered. 


Ill 


TABLE  OF  CONTENTS 


Abstract  iii 

List  of  Illustrations  vii 

1.  INTRODUCTION  1 

2.  DETECTION  AND  COPY  STATISTICS  3 

2.1  Signal  Model  3 

2.2  Maximum  Likelihood  Formulation  4 

3.  REFORMULATION  OF  DETECTION  AND  COPY  7 

3.1  Detection  7 

3.2  Copy  7 

4.  PERFORMANCE  ASSESSMENT  9 

4.1  SNRs  and  Emitter  Separations  9 

4.2  Simulation  Results:  Detection  12 

4.3  Simulation  Results;  Copy  14 

5.  SUMMARY  23 

APPENDIX  A  -  STATISTICAL  SIGNAL  PROCESSING  25 

A.l  Signal  Model  25 

A. 2  Maximum  Likelihood  Decision  Statistic  26 

A. 3  Special  Case:  J  -  I  =  M  29 

APPENDIX  B  -  A  MAXIMIZATION  PROBLEM  31 

APPENDIX  C  -  ARRAY  BEAMSHAPES  33 

APPENDIX  D  -  PERFORMANCE  35 

D.l  Normal  Forms  35 

D.2  Detection  36 

D.3  Special  Case  38 

D.4  Copy  44 


V 


LIST  OF  ILLUSTRATIONS 


Figure 

No. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 


D-1 


Page 


ASNR  loss  in  cochannel  interference.  13 

Detection.  13 

False  alarm  probabilities  (X  >  iV  and  J  =  1  =  M).  15 

Detection  with  q-data  contamination:  p-data  TASIR  30  dB.  15 

Detection  with  q-data  contamination;  p-data  TASIR  20  dB.  16 

Detection  with  q-data  contamination;  p-data  TASIR  20  dB.  16 

Detection  with  q-data  contamination;  p-data  TASIR  10  dB.  17 

Achieved  and  ideal  ASIRs.  17 

Copy  with  q-data  contamination;  p-data  ASIR:  10  dB.  19 

Copy  with  q-data  contamination;  p-data  ASIR:  10  dB.  19 

Copy  with  q-data  contamination;  p-data  ASIR:  20  dB.  20 

Copy  with  q-data  contamination;  p-data  ASIR:  20  dB.  20 

Copy  with  q-data  contamination;  p-data  ASIR:  40  dB.  21 

Copy  with  q-data  contamination  and  0.01  beamwidth  steering  vector  mis¬ 
match;  p-data  ASIR:  10  dB.  21 

Comparison  of  Incomplete  Beta  function,  bounds,  and  asymptotics.  44 


vii 


1.  INTRODUCTION 


Adaptive  detection  and  copy  of  temporally  coherent  (i.e.,  known)  narrowband  waveforms 
using  multiple  sensors  has  been  treated  extensively  in  the  radar  literature.  Kelly  [1]  provides  a 
maximum  likelihood  (ML)  formulation  of  the  problem  (see  also  Kelly  and  Forsythe  [2]).  For  radar 
problems,  a  fundamental  assumption  is  that  the  sensor  response  to  the  signal  of  interest  (SOI) 
is  known.  Variations  on  this  assumption  lead  one  to  study  the  performance  loss  associated  with 
imperfectly  known  sensor  responses  (wavefront  mismatch). 

In  contrast,  the  problem  of  interest  here  is  the  detection  and  copy  of  a  coherent  waveform  with 
a  completely  unknown  wavefront.  Lack  of  knowledge  about  the  wavefront  can  be  due  to  unresolved 
multipath  or,  more  simply,  to  an  uncalibrated  (or  uncalibratable)  antenna  array.  The  problem  of 
wavefront  mismatch  is  absent  given  the  techniques  discussed  in  this  report,  but  for  many  appUca- 
tions  there  can  be  other  problems  associated  with  imperfect  use  of  the  coherent  signal  waveform 
(waveform  mismatch).  The  impact  of  such  problems  is  studied  here.  To  make  the  techniques  more 
concrete,  a  particular  frequency  domain  application  is  discussed.  For  this  application,  the  signal 
processing  is  called  frequency  feature  processing  (FFP). 

In  broad  terms,  the  work  presented  here  specializes  the  results  of  Kelly  and  Forsythe  [2]  to  the 
problems  of  detection  and  copy  of  a  coherent  waveform  with  unknown  wavefront.  Specific  forms  of 
the  detection  statistic  and  adaptive  beamformer  are  presented.  These  forms  are  suited  to  typical 
applications.  Problems  associated  with  waveform  mismatch  are  addressed. 

More  specifically,  Section  2.1  recalls  the  signal  model  of  Kelly  and  Forsythe  [2].  To  simplify 
the  exposition,  the  treatment  of  detection  and  copy  in  the  body  of  this  report  is  expressed  in 
terms  of  a  specific  frequency  domain  application.  However,  the  appendices  provide  a  more  general 
treatment.  Section  2.2  expresses  the  generalized  likelihood  ratio  detection  statistic  and  copy  weight 
of  Kelly  and  Forsythe  [2]  for  the  signal  model  of  interest.  Section  3  reformulates  the  detection 
statistic  and  copy  weight  in  forms  more  appropriate  for  applications.  Section  4  discusses  some 
difficulties  with  detection  and  copy  in  practical  applications  and  quantifies  performance. 
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2.  DETECTION  AND  COPY  STATISTICS 


2.1  Signal  Model 

Let  Z  —  (Z(l), . . . ,  Z{L))  be  an  iV  X  L  array  of  vector  samples  of  the  output  of  an  N  element 
antenna  array.  Each  column  of  Z  represents,  for  FF?  applications,  a  snapshot  of  the  array  output 
in  a  particular  frequency  cell.  It  is  assumed  below  that  these  snapshots  are  independent.  The 
entries  of  Z  are  taken  to  be  jointly  complex  circular  Gaussian.  A  more  general  signal  model  is 
discussed  in  Appendix  A.l. 

One  can  partition  the  data  Z  into  the  N  X  M  subarray  Zp  (p-data)  and  the  N  X  {L  —  M) 
subarray  Z,  (q-data):  Z  =  (Zp  Z,).  It  is  assumed  that  L  >  M  +  N .  Each  column  of  Z  is  assumed 
to  have  the  same  covariance  R.  That  is, 

R  =  E[(Z(fc)  -  E[Z(fc)])  (Z(k)  -  E[Z(fc)])«] 

where  E[-]  denotes  expectation.  However,  the  modeled  mean  of  Z,  is  assumed  to  be  zero  while  that 
of  Zp  is 


VA  Hs 
0  Hns 

under  the  signal  present  (Hs)  and  signal  absent  (Hns)  hypotheses.  F  is  a  A  x  7  array,  the  columns 
of  which  represent  the  array  responses  (signal  wavefronts)  of  J  emitters.  AisaJxM  array  of 
signal  amplitudes. 

The  true  mean  of  Z  need  not  agree  with  the  model.  It  can  be  written  in  general  as 

E[Zp]  =  VAp 
E[Z,]  =  VAg 

where  A  =  (ApAg)  is  a  7  x  T  array  expressing  contributions  of  the  signal  to  both  the  p-data  and 
q-data.  V  is  a  A  x  7  array  expressing  the  true  array  responses  of  7  emitters.  When  it  is  important 
to  make  a  distinction,  a  script  font  is  used  to  signify  true  properties  of  the  signal  mean  as  opposed 
to  properties  of  the  model  (expressed  in  a  roman  font).  Signal  presence  in  the  q-data  is  a  form  of 
model  mismatch  that  impacts  the  detection  and  copy  procedures  described  below. 

For  FFP  applications  (the  treatment  is  more  general  in  the  appendices)  7  =  1.  Then  A  is 
a  vector  representing  the  signal  components  in  the  p-data  and  q-data.  Signal  in  the  q-data  is  a 
type  of  modeling  inaccuracy  due,  for  example,  to  uncertainty  in  the  signal  bandwidth.  For  FFP 
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detection,  Af  =  1  is  assumed  in  addition.  Then  the  modeled  .4  is  a  scalar  that  lives  in  one  and  only 
one  frequency  cell,  which  may  or  may  not  be  known  a  priori. 

The  model  of  V  may  vary  considerably,  leading  to  different  forms  of  parameter  estimators. 
V  is  assumed  to  be  completely  unknown  in  the  FFP  applications.  Thus  the  procedures  described 
below  can  be  used  with  uncalibrated  arrays.  Note  that  there  is  no  modeling  mismatch  associated 
with  V  in  this  case;  all  V’s  are  included  in  the  model. 

An  example  of  a  waveform  amenable  to  FFP  is  a  narrowband  AM  signal  with  significant 
energy  in  the  carrier.  For  detection  purposes,  the  signal  is  modeled  as  a  pure  tone  at  the  carrier 
frequency.  Thus  detection  suffers  from  model  inaccuracies  due  to  the  true,  multicell  signal  band¬ 
width.  The  effect  of  this  modeling  error  is  assessed  in  Section  4.2.  Copy  also  relies  on  the  tone 
model,  but  a  guard  band  can  be  introduced  to  address  the  effects  of  an  inaccurate  model  (see 
Section  4.3). 

2.2  Maximum  Likelihood  Formulation 

Letting  5  ,  the  generalized  likelihood  ratio  (see  Appendix  A. 2  for  a  concise  descrip¬ 

tion  of  the  test  in  a  more  general  situation  as  well  as  in  the  special  case  treated  here)  for  deciding 
between  the  hypotheses  Hs  and  Hns  is 


v^s-^v 

V»iZZ^)-W 

The  associated  copy  weight  is  expressed  by 


(1) 


(2) 


The  expression  W^Zp  provides  a  ML  estimate  of  the  signal  amplitude  A. 

To  complete  the  detection  statistic  and  determine  the  copy  weight  W,  Equation  (1)  must  be 
maximized  over  V.  V  can  be  modeled,  for  example,  by  a  calibration  table  for  the  antenna  array. 
Because  detection  and  copy  with  uncalibrated  arrays  are  of  principal  interest  here,  V  is  assumed 
to  be  completely  unknown.  Thus  Equation  (1)  must  be  maximized  over  all  nonzero  V .  To  do  this, 
one  can  use  the  following  well-known  result,  valid  for  any  hermitian  matrices  A  and  B,  with  B 
positive  definite. 


X^AX  X^B-^I'^AB-^I'^X  .  ,0-1/2. 0-1/2,  ,  /o-i., 

BX  ~  AB  )  —  Anaax(-®  ^4), 

where  the  maximization  occurs  over  vectors  X  and  Ai„ax(  )  denotes  the  largest  eigenvalue  of  its 
argument.  Let  enjax(-)  denote  the  principal  eigenvector  of  its  argument.  Then 
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-X^max  =  €nvax(5 

Applying  the  above  fact  to  Equation  (1),  one  has 
•wH  c— 1  y 

mw  =  A^[(5  +  2,Z,")S-'l  =  1  +  Z«5-'Zp.  (3) 

which  occurs  when 

14iax  =  CnvaxK'S’  +  ZpZ^)S  ']  =  Zp. 

The  copy  weight  associated  with  Vniax  is 

W  =  S-^Zp{Z^S-^Zp)-\  (4) 

Detection  based  on  Equation  (3)  enjoys  a  strong  constant  false  alarm  rate  (CFAR)  property. 
Specifically,  the  distribution  of  the  detection  statistic  under  does  not  depend  on  the  spatial 
distribution  of  interference. 
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3.  REFORMULATION  OF  DETECTION  AND  COPY 


3.1  Detection 

The  detection  and  copy  statistics  presented  above  have  more  useful  forms,  especially  when 
the  frequency  cell  containing  the  SOI  is  not  known  but  must  be  located  as  part  of  the  detection 
procedure.  The  new  statistics  are  derived  by  applying  the  Sherman-Morrison-Woodbury  identity 
(see  Golub  and  Van  Loan  [3]) 

(5+  ZpZp*)-’  =  5-^  -  5-'Zp(l  +  S-^ Z» S~y 

Consider  detection  first.  The  detection  statistic  is  given  by  Z^S~^Zp.  This  applies  to  a 
specific  frequency  cell.  If  the  cell  containing  the  SOI  is  not  known,  the  covariance  estimate  (unnor¬ 
malized)  5  must  be  reformulated  for  each  cell  tested.  However, 


Zp"(ZZ")-'Zp  =  Zp"(5-h  ZpZp")-'Zp  = 


Z»S-^Zp 
l  +  Z»S-^Zp' 


Because  a:/(l  +  x)  is  a  monotonic  increasing  function,  an  equivalent  form  of  the  detection  statistic 
is  given  by  the  left-hand-side  of  the  above  equation,  namely 

Zp«(ZZ«)-'Zp.  (5) 

This  version  of  the  statistic  uses  only  one  unnormalized  covariance  estimate  ZZ^  for  all  cells  tested. 

3.2  Copy 

A  similar  (though  less  useful)  reformulation  is  possible  for  the  copy  weight  baised  on  a  partic¬ 
ular  cell: 

(ZZ")-^Zp(Z"(ZZ")-^Zp)-'  =  [5-‘Zp(l-hZp^5-^Zp)-^] 

•[(Zp^5-'Zp)(l  -h  Zp"5-'Zp)-']-^ 

=  5-'Zp(Zp"S-^Zp)-^ 

The  right-hand  side  of  the  above  equation  can  be  recognized  as  the  adaptive  copy  weight  introduced 
in  Equation  (4).  Thus  the  adaptive  weight  can  be  written  in  the  form 


(ZZ^)-iZp(Zf(ZZ^)-'Zp)-'. 


(6) 
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4.  PERFORMANCE  ASSESSMENT 


Both  the  detection  statistic  of  Equation  (5)  and  copy  weight  of  Equation  (6)  can  suffer  from 
modeling  inaccuracies  in  practical  applications.  The  signal  model  assumes  that  the  SOI  is  concen¬ 
trated  in  a  single  frequency  cell  while  the  interference,  in  principle,  is  found  in  all  cells.  Typical 
applications  of  FFP  involve  signals  with  spectral  lines  due  to  residual  carriers  or  due  to  periodicities 
in  the  data.  Thus  a  substantial  fraction  of  the  SOI  power  is  concentrated  in  one  or  several  cells,  but 
there  can  be  considerable  signal  power  in  other  cells  as  well.  In  order  to  characterize  performance, 
these  modeling  inaccuracies  must  be  taken  into  account. 

4.1  SNRs  and  Emitter  Separations 

Typically,  the  performance  of  adaptive  array  algorithms  is  assessed  in  terms  of  emitter  levels 
and  emitter  separations.  Some  basic  definitions  of  both  are  introduced  next.  More  discussion  of 
emitter  separations  is  presented  in  Appendix  C. 

The  array  signal-to-noise  ratio  (ASNR)  of  an  emitter  in  a  background  environment  (noise  and 
possibly  some,  but  not  all,  interference)  is  the  SNR  at  the  output  of  an  optimal  beamformer.  It  is 
the  largest  SNR  achievable  with  adaptive  beamforming  in  the  given  environment.  What  constitutes 
the  “background  environment”  as  opposed  to  “interference”  is  largely  a  question  of  what  is  known 
as  opposed  to  unknown.  Because  the  background  is  often  known,  it  is  incorporated  i’.  the  definition 
of  ASNR  because,  in  this  case,  the  data  can  be  prewhitened  by  the  background  before  adaptive 
processing,  which  handles  additional,  unknown  interference.  In  other  words,  it  is  useful  to  include 
any  known  background  in  the  definition  of  ASNR  in  order  tc  gauge  performance  against  unknown 
(typically  discrete)  sources.  In  some  cases  it  is  convenient  to  exclude  known  interference  from  the 
concept  of  background  environment.  For  example,  the  noise  floor  may  be  determined  by  spatially 
diffuse  environmental  sources  rather  than  by  thermal  noise.  Such  a  background,  if  it  does  not 
appear  spatially  white,  can  be  suppressed  to  some  extent  by  adaptive  techniques.  By  excluding 
such  interference  from  the  definition  of  background,  this  suppression  can  be  assessed  [in  terms  of 
array  signal-to-interference  ratio  (ASIR);  see  below].  However,  it  is  worth  noting  that  adaptive 
techniques  provide  the  most  benefit  against  discrete  sources. 

The  calculation  of  ASNR  is  based  on  the  covariance  decomposition 
signal  cov. 

R=  Rs  +  Rn 

background  cov. 


In  the  following,  the  reader  can  safely  assume  Rn  =  In->  where  In  is  the  N  x  N  identity  matrix, 
representing  a  spatially  white  background.  This  leads  to  the  most  common  definition  of  ASNR;  but 
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as  the  discussion  above  indicates,  other  definitions  can  be  useful.  Adaptive  beamforming  maximizes 
the  SNR  at  the  output  of  the  beamformer,  forming 


W^RsW  . 

W^^RnW  ~  ^ 


which  is  achieved  when 

^max  ~  Rs^- 

Note  that  there  is  assumed  to  be  no  correlation  between  signal  and  background.  In  practice, 
Rs  =  ■,  expressing  the  contribution  of  a  single  coherent  emitter  to  the  total  covariance.  V 

(a  column  vector  here)  is  the  emitter  steering  vector  and  p5  is  the  emitter  power.  In  this  case, 

ASNR  =  \r^{R-f,^Rs)  =  PS^^^{R~J,^'^VV»R-^^''^)  =  psiV^R^^V) 


and 


v^max  =  R~nV. 

In  an  interference  environment,  the  covariance  is  expressed  by  i?  =  Rs-\-Ri,  where  Ri  >  Rn  is 
the  contribution  of  interference  sources  (typically  discrete)  as  well  as  the  background  environment. 
With  the  signal  model  discussed  above,  the  optimal  ASIR  is  given  by  psiV^ RJ^V).  This  leads 
one  to  define  the  interference  induced  loss 

V^RjW 

V^RJf^V 

This  is  the  loss  of  signal-to-interference  ratio  (SIR)  due  to  the  presence  of  the  additional  interference 
contained  in  R/. 

In  a  special  case  of  interest,  the  added  interference  is  due  to  a  single  coherent  emitter,  so 
that  Ri  =  Rs  +  pfWW^ ,  with  W  the  steering  vector  of  the  interferer.  Define  the  separation  in 
beamwidths  between  the  two  steering  vectors  V  and  W  as 


b  =  (2/7r)arcos( 


{v^R-s'wi 


(7) 
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where  {V^ etc.  Appendix  C  explains  this  definition.  When  Rn  k  /at,  this 

definition  agrees  closely  with  the  physical  notion  of  beamwidth  for  values  of  b  less  than  one-half. 
Note  that  6=1  corresponds  to  the  orthogonality  of  vectors  V  and  W.  Because  (Sherman-Morrison- 
Woodbury  again;  see,  for  example,  Golub  and  Van  Loan  [3]) 


=  V 


PI 


l+pjWffRj/W 


R]yWW 


H 


Rj^ , 


the  interference  induced  loss  becomes 


piw^Riyw  \v^Rjywy 

l-]-piW»R]yW  ■  {W^iR]yW)iV»R]yV) 


ASNR/ 

1  -)-  ASNR/ 


cos^(7r6/2), 


(8) 

(9) 


where  ASNR/  is  the  ASNR  of  the  interferer.  As  ASNR/  increases,  this  expression  reaches  the 
asymptote  sm^{nb/2).  Figure  1  shows  some  examples  of  the  interference  induced  loss  as  a  function 
of  ASNR/  and  the  beamwidth  separation  of  the  emitters.  Note  that  an  emitter  separation  of 
one-half  beamwidth  leads  to  a  loss  (the  reciprocal  of  Equation  (9)  is  plotted  on  dB  scale,  for 
convenience)  of  at  most  3  dB.  If,  on  the  other  hand,  the  emitters  are  within  0.1  beamwidths  of  each 
other,  the  loss  can  approach  16  dB. 

Some  additional  comments  on  SIRs  are  necessary  for  the  discussion  below.  The  definition  of 
ideal  ASIR  (when  J  =  1,  M  arbitrary)  can  be  stated  in  terms  of  the  signal  model  of  Section  2.1 
using  the  terminology  appropriate  for  FFP  applications.  A  more  general  definition  is  given  in  the 
appendices. 

Let  Ak,  the  fcth  component  of  A,  be  the  signal  amplitude  in  the  A:th  frequency  cell,  then  the 
ideal  ASIR  becomes  ||A]|^  (V^ R~^V)/L.  Note  that  this  is  really  an  average  ideal  ASIR  taken  over 
L  observations.  This  can  also  be  called  the  p-data  ideal  ASIR  when  ^4*;  =  0  in  the  q-data.  When 
there  is  SOI  contamination  of  the  q-data,  one  can  define  separate  p-data  and  q-data  ideal  ASIRs 
by  restricting  ||>41P  to  the  p  or  q  components  (|!>Ip|p  or  |l^g|P). 

There  is  a  distinction  between  the  ideal  ASIR  just  defined,  which  is  a  per  sample  SIR,  ap¬ 
propriate  for  characterizing  copy  performance,  and  the  ideal  total  ASIR  (TASIR),  which  is  needed 
to  characterize  detection.  The  ideal  TASIR,  for  FFP  detection  applications,  is  L  times  the  ideal 
ASIR.  Recall  that  L  is  the  number  of  available  samples.  The  ideal  TASIR  incorporates  coherent 
and  incoherent  (M  >  1)  SOI  integration  into  the  SIR.  For  convenience,  the  ideal  TASIR  is  called 
simply  TASIR  below. 
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4.2  Simulation  Results:  Detection 

All  the  plots  in  this  section  are  based  on  the  TASIR  (or  sometimes  the  p-data  and  q-data 
versions). 

Figure  2  shows  the  cumulative  distribution  functions  (CDF)  of  the  normalized  detection  statis¬ 
tic  Equation  (10)  for  SOI  TASIRs  of  10,  20,  and  30  dB.  The  figure  is  based  on  a  Monte  Carlo 
simulation  utilizing  the  statistical  normal  forms  discussed  in  Appendix  D.l.  The  dimensional  pa¬ 
rameters  are  IV  =  4,  M  =  1,  and  L  =  1000.  Also  shown  in  the  figure  is  the  statistic  with  no  signal 
present  (false  alarm  statistic).  Note  that  the  detection  statistics  axe  essentially  constant  when  the 
SOI  TASIR  is  much  larger  than  unity. 

To  see  how  the  figures  can  be  used,  consider  a  signal  with  an  ASNR  of  -4  dB  that  is  0.1 
beamwidth  away  from  a  very  strong  interferer.  Figure  1  shows  that  the  interference  induced  loss 
is  at  most  16  dB.  Thus  the  ideal  ASIR  is  at  least  —20  dB.  With  L  =  1000,  the  TASIR  becomes  10 
dB.  The  corresponding  curve  in  Figure  2  can  be  used  to  determine  detection  performance.  Note 
that  the  signal  can  be  detectable  even  if  the  ASIR  does  not  support  good  copy. 

The  false  alarm  probability  can  be  characterized  when  J  =  1  =  M.  See  Appendix  D.3.2 
for  the  calculations  and  some  simple  approximations.  Figure  3  shows  false  alarm  probability  for 
=  2, 4, 8, 16  when  L  is  large.  The  normalized  detection  statistic  of  Equation  (10)  is  used.  Because 
Figure  3  shows  a  per  look  false  alarm  probability,  very  low  values  are  displayed.  The  false  alarm 
rate  is  independent  of  the  spatial  distribution  of  the  interference.  However,  it  does  depend  on  the 
spectral  distribution  of  the  interference.  If  the  interference  is  not  spectrally  flat,  higher  false  alarm 
rates  can  result. 

Appendix  D.3  provides  an  approximation  for  the  normalized  detection  statistic 


L~M 

M 


[z^s-^z,] 


(10) 


when  L  1.  Assume  M  =  1,  because  this  is  the  case  of  most  interest  for  FFP  applications.  Let 
TASiRp  denote  the  TASIR  of  the  SOI  in  the  p-data.  Similarly,  let  TASIR,  denote  the  TASIR  of  the 
SOI  in  the  q-data.  The  detection  statistic  is  essentially  constant  when  TASIRp  1  and  assumes 
the  approximate  value 


TASiRp 
-I- TASIR,' 


(11) 


When  TASIR,  C  i,  the  approximate  value  is  TASIRp.  When  TASIR,  >  L,  this  value 
saturates  at  about  TTASIRp/TASIR,.  In  either  case,  a  significant  peak  in  the  statistic  occurs  at 
lines  in  the  SOI  spectrum  because,  in  this  case,  the  value  of  TASIRp/TASIR,  is  comparable  with 
unity  and  because  L  is  large. 
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Figure  4  is  based  on  the  same  dimensional  parameters  as  Figure  2.  The  SOI  has  a  fixed  (for 
each  curve)  TASIR  of  30  dB  in  the  p-data  and  TASIRs  of  30,  20,  and  10  dB  in  the  q-data.  Note 
that  the  SOI  contamination  of  the  q-data  leads  about  a  3  dB  loss  in  the  detection  statistic  given  30 
dB  contamination,  consistent  with  the  approximations  above.  Contamination  of  20  or  10  dB  leads 
to  little  loss  in  detectability. 

Similar  comments  apply  to  Figure  5  and  Figure  6  where  the  p-data  TASIR  is  reduced  to 
20  dB,  the  q-data  TASIRs  become  20,  10,  and  0  dB,  and  the  number  of  looks  is  1000  (Figure  5) 
and  100  (Figure  6). 

Finally,  Figure  7  shows  an  SOI  with  10  dB  p-data  TASIR  and  q-data  TASIRs  of  10,  0,  and 
— 10  dB.  There  are  1000  looks.  Detectability  suffers  no  essential  loss. 

4.3  Simulation  Results:  Copy 

The  copy  performance  achieved  with  the  adaptive  beamformer  of  Equation  (4)  is  expressed 
in  terms  of  the  ASIR  at  the  output  of  the  copy  beam.  In  Appendix  D  it  is  shown  that  this  achieved 
ASIR  (which  is  a  random  variable)  has  a  distribution  that  depends  only  on  the  ideal  ASIR  and 
on  the  dimensional  parameters  X  and  N  (recall,  Af  =  1  in  the  applications  discussed  here).  In 
other  words,  the  achieved  SIR  depends  on  a  particular  scenario  through  the  ideal  ASIR.  Scenario 
details  such  as  interference  spatial/polarization  distributions  are  irrelevant  to  characterizing  copy 
performance  unless  they  alter  the  ideal  ASIR. 

Figure  8  shows  the  ideal  and  achieved  ASIRs  for  a  four  element  array  using  100  samples. 
The  vertical  line  in  the  figure  indicates  ideal  ASIR,  namely  10  dB.  This  line  also  represents  the 
CDF  of  an  ideal  adaptive  beamformer.  The  curved  line  shpws  the  CDF  of  achieved  ASIR  (in  dB) 
provided  by  FFP.  The  achieved  ASIR  is  within  a  few  tenths  of  a  dB  of  the  ideal.  This  is  due  to  the 
large  number  of  samples  (100  here)  typical  of  an  FFP  application  with  a  consequent  large  coherent 
integration  gain.  In  this  example,  the  ideal  TASIR  is  30  dB. 

When  signal  is  present  in  the  q-data,  performance  is  degraded.  As  shown  in  Appendix  D.4, 
the  achieved  ASIR  depends  only  on  the  ideal  ASIRs  of  the  SOI  in  both  the  p-data  and  q-data, 
as  well  as  on  the  dimensional  parameters  mentioned  above.  The  next  set  of  figures  shows  the 
performance  degradation  as  a  function  of  these  ideal  ASIRs.  Each  figure  is  based  on  a  fixed  ideal 
ASIR  (shown  in  the  caption)  for  the  SOI  in  the  p-data.  The  ideal  ASIR  of  the  SOI  in  the  q-data  is 
varied  as  indicated  in  the  legend.  The  term  “No  Q-DATA”  means  that  the  signal  lies  only  in  the 
p-data.  The  dimensional  parameters  are  varied  as  shown  on  the  figures. 
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CUMULATIVE  DISTRIBUTION  FALSE  ALARM  PROBABILITY 
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Figure  3.  False  alarm  probabilities  (L'^  N  and  J  ~  \  =■  M ). 
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Figure  4-  Detection  with  q-data  contamination;  p-data  TASIR  SO  dB. 
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CUMULATIVE  DISTRIBUTION  CUMULATIVE  DISTRIBUTION 
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Figure  5.  Detection  with  q~data  contamination;  p-data  TASIR  20  dB. 
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Figure  6.  Detection  with  q-data  contamination;  p~data  TASIR  20  dB. 
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FRACTION  CUMULATIVE  DISTRIBUTION 
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Figure  7.  Detection  with  q-data  contamination;  p-data  TASIR  10  dB. 
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Figure  8.  Achieved  and  ideal  ASIRs. 


17 


The  reader  should  note  tb^t  the  CDFs  of  achieved  ASIR  shown  below  are  based  on  the  p-data 
ideal  ASIR.  This  is  a  convenient  way  to  present  the  results  given  some  of  the  processing  options 
described  below  (i.e.,  guard  bands).  However,  copy  performance  involving  the  p-data  is  not  directly 
of  interest  because  that  component  of  the  signal  is  already  known  upon  detection  and,  hence,  does 
not  contain  any  additional  information.  It  is  more  likely,  but  not  necessary  (see  below),  that  the 
achieved  ASIR  of  the  q-data  is  of  interest.  In  this  case,  the  CDF  curves  should  be  translated  by 
the  difference  (in  dB)  between  the  q-data  and  p-data  ideal  ASIRs.  For  example,  a  CDF  curve  of 
achieved  ASIR  based  on  an  ideal  p-data  ASIR  of  30  dB  and  an  ideal  q-data  ASIR  of  20  dB  should 
be  translated  left  by  10  dB. 

Figure  9  shows  that  there  is  essentially  no  loss  in  achieved  ASIR  when  X  >  1  unless  the 
q-data  ASIR  is  much  greater  than  the  p-data  ASIR.  Here,  the  ideal  p-data  ASIR  is  low  (10  dB). 
The  loss  due  to  q-data  contamination  increases  as  the  number  of  looks  (X)  is  reduced  (Figure  10) 
and  as  the  p-data  ideal  ASIR  is  increased  (Figure  11  and  Figure  12).  However,  these  losses  are  not 
significant  as  long  as  the  q-data  ASIR  is  much  less  than  the  p-data  ASIR. 

If,  as  mentioned  above,  the  achieved  ASIR  for  the  q-data  is  of  interest,  one  can  rephrase  the 
preceding  statements.  As  long  as  the  q-data  ASIR  is  less  than  (approximately)  the  p-data  ASIR, 
copy  performance  improves  with  increasing  q-data  ASIR  (but  fixed  p-data  ASIR).  Once  the  q-data 
ASIR  exceeds  the  p-data  ASIR,  copy  performance  degrades.  For  example,  Figure  11  shows  that  a 
q-data  ASIR  of  20  dB  leads  to  a  median  achieved  ASIR  (for  the  q-data)  of  about  18  dB.  However, 
a  q-data  ASIR  of  30  dB  (which  exceeds  the  p-data  ASIR  by  10  dB)  leads  to  a  median  achieved 
ASIR  of  about  15  dB.  In  other  words,  copy  performance  has  been  degraded. 

It  is  worth  bearing  in  mind  that  the  copy  performance  described  in  the  figures  is  conditioned 
on  detecting  the  SOI.  When  this  is  not  possible  with  high  probability,  copy  performance  can  deviate 
significantly  from  that  shown  in  the  figures  due  to  mislocating  the  cell  containing  the  SOI  spectral 
line  (for  FFP  applications).  In  the  figures  shown  here  only  the  30  dB  q-data  curve  of  Figure  12  is 
effected. 

It  is  possible  for  the  SOI  contamination  of  the  q-data  to  have  a  more  serious  impact  on 
performance.  Figure  13  shows  the  effect  of  contamination  at  a  high  ideal  ASIR.  Figure  14  examines 
the  same  scenario  with  the  addition  of  a  0.01  beamwidth  mismatch  between  the  SOI  wavefront  in  the 
q-data  and  the  p-data  (which  is  taken  to  be  truth).  Note  that  the  losses  due  to  SOI  contamination 
are  considerably  larger.  Steering  vector  mismatch  can  arise  for  a  number  of  reasons.  For  example, 
large  multipath  delays  can  cause  wavefronts  to  vary  across  the  signal  bandwidth,  violating  the 
assumption  that  the  SOI  is  narrowband. 
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FRACTION  FRACTION 
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Figure  9.  Copy  with  q-daia  contamination;  p~data  ASIR:  10  dB. 
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Figure  10.  Copy  with  q-data  contamination;  p-data  ASIR:  10  dB. 
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FRACTION  FRACTION 
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Figure  13.  Copy  with  q-data  contamination;  p-data  ASIR:  40  dB. 
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Figure  14-  Copy  with  q-data  contamination  and  0.01  beamwidth  steering  vector  mis¬ 
match;  p-data  ASIR:  10  dB. 
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For  FFP  applications,  the  q-data  consists  of  all  frequency  cells  ex  Jt  the  cell(s)  containing 
the  SOI  lines.  Thus  the  ideal  ASIR  in  the  q-data  can  be  large.  This  could  severely  degrade  copy 
performance,  particularly  when  there  is  some  wavefront  variation  over  the  band,  as  the  figures 
indicate.  To  avoid  this  problem,  a  guard  band  can  be  centered  about  the  line(s).  One  can  think  of 
preprocessing  the  data  by  notching  out  everything  in  the  guard  band,  but  allowing  the  line(s)  to 
come  through.  The  preprocessed  data  does  not  have  significant  SOI  contamination  in  its  q-data  if 
the  guard  band  is  chosen  appro  -ately.  Of  course,  SOI  power  is  thrown  away  (but  only  for  weight 
determination)  if  the  guard  band  is  notched  out.  Instead,  one  can  redefine  the  p-data  so  that  it 
contains  the  guard  band  instead  of  just  the  iine(s).  However,  if  the  line(s)  contains  a  large  fraction 
of  the  signal  power,  this  will  not  improve  performance  significantly.  Of  course,  if  a  guard  band 
is  used,  the  calculation  of  the  q-data  (instead  of  p-data)  achieved  ASIR  described  above  must  be 
modified. 
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5.  SUMMARY 


The  work  in  Kelly  and  Forsythe  [2]  is  reformulated  and  extended  for  applications  involving 
the  detection  and  copy  of  coherent  waveforms  with  unknown  wavefronts.  An  important  application 
of  this  work  (and  a  convenient  example  in  the  exposition)  is  called  frequency  feature  processing 
(FFP).  For  this  application,  the  waveforms  are  narrowband  with  significant  power  concentrated  in 
one  or  several  spectral  lines.  These  waveforms  have  unknown  wavefronts  that  differ  from  those  of 
any  cochannel  interferers.  The  wavefront  differences  allow  spatial  nuDing  of  the  interference. 

Detection  and  copy  performance  are  discussed.  Direction-finding  is  not  considered  because 
antenna  calibration  is  assumed  to  be  unavailable. 

Detection  enjoys  a  constant  false  alarm  rate  property  such  that  the  probability  of  false  alarm 
does  not  depend  on  the  spatial  distribution  of  the  environment,  including  interference.  It  must,  of 
course,  depend  on  the  spectral  distribution  of  the  environment /interference.  Detection  performance 
relies  on  spatial  nulling  in  addition  to  coherent  integration  of  the  signal  of  interest’s  (SOI)  spectral 
line(s).  As  a  result,  signal  detectability  is  substantially  greater  than  that  achievable  with  a  discrete 
Fourier  transform  (DFT)  alone 

A  typical  application  involves  a  SOI  with  spectral  line(s)  (i.e.,  carrier  or  rat"?  lines)  in  cochan¬ 
nel  broadband  (spectrally  flat)  interference.  Because  the  SOI  can  have  significant  power  outside  of 
the  lines,  SOI  self-nulling  is  a  potential  problem.  Self-nulling  can  occur  when  the  SOI  signal  level 
is  large  and  most  of  the  SOI  power  is  not  contained  in  the  line(s).  Loss  of  performance  due  to  this 
problem  is  quantified. 

The  appendices  discuss  a  more  general  signal  model  than  that  treated  in  the  body  of  the 
report.  Applications  of  this  more  general  model  are  not  pursued  here. 
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APPENDIX  A 

STATISTICAL  SIGNAL  PROCESSING 


The  treatment  in  this  appendix  is  more  general  than  that  in  the  body  of  the  report.  In 
particular,  the  signal  model  is  much  more  elaborate.  After  recalling  the  maximum  likelihood 
statistic  in  a  form  given  in  Kelly  and  Forsythe  [2],  some  interesting  reformulations  are  presented. 

A.l  Signal  Model 

Let  Z  he  dm  N  y.  L  array  of  vector  samples  of  the  output  of  an  N  element  antenna  array. 
Each  column  of  Z  is  a  snapshot  of  the  array  output  at  a  particular  time.  It  is  assumed  below  that 
these  snapshots  are  independent  although  the  initial  formulation  is  somewhat  more  general.  The 
entries  of  Z  are  taken  to  be  jointly  complex  circular  Gaussian  with  mean 

E[Z]  =  VAT.  (A.l) 

V  isa.  N  X  J  array,  the  columns  of  which  represent  the  array  responses  of  J  emitters.  T  is  &  M  x  L 
array,  the  rows  of  which  can  be  interpreted  as  time  domain  basis  waveforms.  A  is  a.  J  x  M  array 
of  signal  amplitudes.  In  one  case  of  interest  here  M  =  1  and  J  =  1.  Then  VAT  expresses  the  time 
history  of  a  single  coherent  waveform.  In  the  formulation  of  hypothesis  testing  problems  in  Kelly 
and  Forsythe  [2],  the  A  array  is  assumed  to  be  completely  unknown.  However,  prior  knowledge  of 

V  and  T  may  vary  considerably,  leading  to  different  forms  of  parameter  estimators. 

The  covariance  of  Z  is  given  by 

(cov(Z))(,, E[{Zii  -  Wi)  {Zki  -  Wi)*]  (A.2) 

where  the  overbar  designates  expected  value  and  the  *  denotes  complex  conjugate.  It  is  often  the 
case  that  the  covariance  has  the  special  structure  expressed  by 

cov(Z)  =  il  0  C*  (A.3) 

where  the  Kronecker  (tensor)  product  of  matrices  A  and  B  is  defined  by 

(A  0  —  ^ikSjl  (A.4) 
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The  tensor  product  structure  of  the  covariance  can  be  interpreted  as  a  factoring  into  spatial 
(i?)  and  temporal  (C)  factors.  For  the  applications  of  interest  here,  the  factor  C  is  assumed  to  be 
known.  It  can  express,  for  example,  the  time  domain  correlation  imposed  on  a  wideband  waveform 
by  the  bandpass  characteristic  of  a  receiver. 

A.2  Maximum  Likelihood  Decision  Statistic 

Some  additional  notation  and  machinery  is  required.  Consider  data-independent  matrices  F 
and  D  that  pre-  and  postmultiply  the  data,  Z.  F  can  represent  spatial  processing  of  the  data  while 
D  can  represent  temporal  processing.  It  is  easy  to  show  that 

^[FZD]  =  F^[Z\D  (A.5) 

and  that  (see  Kelly  and  Forsythe  [2],  Equation  Al'44) 

cow{FZD)  =  {FRF^)  ®  {D^CDf.  (A.6) 

In  the  following,  it  is  assumed  that  C  =  /£,,  the  Lx  L  identity  matrix. 

An  orthonormal  basis  of  the  row  space  of  T  can  be  chosen  by  defining 

P  =  (A.7) 

This  is  an  M  X  X  array,  the  rows  of  which  span  the  same  subspace  as  the  rows  of  T.  Furthermore, 
the  rows  of  P  are  orthonormal,  that  is 

PP^  =  Jm-  (A.8) 

Choose  3iK  =  L  —  MhyL  array,  Q,  the  rows  of  which  are  orthonormal  and  span  the  orthocom¬ 
plement  of  the  row  space  of  P.  The  P  and  Q  arrays  have  the  following  properties: 


/l 

=  p^pfq^q 

PQ^ 

=  0 

ppH 

=  Im 

QQ^ 

=  Ik- 
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Postmultiply  the  data  array  Z  with  a  data  independent,  unitary,  Lx  L  matrix,  Ulj ,  where 


Ul  = 


(A.9) 


Then 


ZUjj  =  (ZP"  ZQ^)  =  (Zp  Z,)  (A.IO) 

is  a  left-right  partition  of  the  transformed  data.  Postmultiplication  by  Ui  constitutes  time-domain 
processing  of  the  data  array,  Z,  with  a  bank  of  orthogonal  filters.  The  covariance  machinery 
presented  above  tells  us  that 

cov(ZUlj)  =  R®  UlIlUH  =  (A.ll) 

In  order  to  find  the  transformed  mean  of  Z,  write 

VAT  =  VBP  (A.12) 

where 

B  (A.13) 


Then 


E[Zp]  =  VBPP^  =  VB  ;  E[Z,]  =  VBPQ^  =  0.  (A.14) 

Thus  Zp  and  Z,  are  independent  Gaussian  arrays  with  independent  columns  and  identical  col¬ 
umn  (spatial)  covariances.  Z,  has  zero  mean  while  Zp  has  a  nonzero  mean  containing  the  signal 
parameters.  Let  S  below.  In  addition,  note  that  ZZ^  —  +  S. 

From  Kelly  and  Forsythe  [2]  (Equations  2-1,  2-2,  2-8,  2-57),  the  generalized  likelihood  ratio 
(GLR)  decision  statistic,  maximized  over  A,  can  be  written 
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[\V»{ZZ«)-W\)  maxRir-NLiii\-L^-trR-^zz«  ’ 

with  the  maximum  likelihood  (ML)  estimate  of  A  given  by  (Kelly  and  Forsythe  [2],  Equation  2-49) 

W^ZT^{TT^)~^ 


where 

w  =  s~'^v{v^s-'^v)-\ 


For  the  special  case  treated  in  the  body  of  the  report, 

T  =  {Im 

and  hence  Z  =  (Zp  Zq),  so  that 

ZT^{TT»r^  =  {Zp  Zq){lM  0M,L-Mf  =  Zp. 

Thus  the  ML  estimate  of  A  is  W^Zp. 

Let  Afc(-)  denote  the  fcth  largest  eigenvalue  of  its  argument.  In  Kelly  and  Forsythe  [2]  Equation 
2-61  it  is  shown  that  (see  also  Appendix  B) 

=  n(l  +  Ai(Z»S-'Z,)).  (A.16) 

fc=l 

The  last  step  follows  from  the  fact  that  the  nonzero  characteristic  values  (including  multiplicity, 
i.e.,  counting  each  occurrence  of  repeated  values)  of  the  products  of  any  nxl  F  and  Ixn  D,  namely 
FD  and  DF,  are  the  same.  To  see  this,  recall  that  the  polynomial  \In  +  xC\,  for  any  square  matrix 
C,  has  degree  rank  C  and  has  zeros  of  the  form  — A~^  where  A  runs  over  the  nonzero  characteristic 
values  (roots  of  the  characteristic  polynomial)  of  C.  In  particular,  if  C  is  diagonalizable,  these 
characteristic  values  are  the  same  as  eigenvalues.  From  a  determinant  identity  given  in  Appendix 
1  (Equation  Al-3)  of  Kelly  and  Forsythe  [2],  one  has  |/„  -|-  xFD\  =  |/i  -|-  xDF\,  and  hence  FD  and 
DF  have  the  same  nonzero  characteristic  values.  In  the  application  above,  both  FD  and  DF  are 
diagonalizable  and  thus  have  the  same  nonzero  eigenvalues. 
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The  result  above  can  be  recast  in  a  more  useful  form  for  applications.  Using  the  Sherman- 
Morrison- Woodbury  identity  (see  Golub  and  Van  Loan  [3])  one  can  write 

Z»(ZZ^r^Zp  =  + 

=  Z"(5-'  -  S-'^Z^ilM  +  Zf  5-iZp)-'Zf  5-')Zp 

=  iz^s-^z^){iM  +  z»s-^Zp)-\ 

Because  /(x)  x/(l  +  x)  is  monotonic  increasing  in  x  for  x  >  0,  it  follows  that 


Afc(Z"(ZZ^r'Zp) 


Afc(Zf5-’Zp) 

1  +  Xk(Z»S-^Zp) 


and  hence  that 


1  +  A;,(Z"5-^Z^)  (1  -  A,(Zf  (ZZ"r>Zp))-^ 


(A.17) 


(A.18) 


n(l  +  Ai(Z"S-iZ,)) 

A;=l 

n(l  -  Xk{Z^{ZZ»r^Z^))-\  (A.19) 

fc=i 


\Jm  -  Z"(ZZ")-'Zpr^  (A.20) 

A.3  Special  Case:  J  =  1  =  M 

With  Pz  ^  Z^{ZZ^)~^Z  and  J  =  1,  the  detection  test  statistic  Equation  {A.19)  becomes 
(1  -  \r.^{Z«{ZZ^)-^Z^)r^  =  (1  -  Xn^iPPzP^))-^.  (A.21) 


Thus 


\V”S-^V\ 


When  J  =  M  this  becomes 
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With  M  =  1  this  is  just 


(i-iippzfr^ 

Thus,  if  P  varies  over  a  parameterized  family  of  waveforms,  the  GLR  detection  test  involves  max¬ 
imizing 

\\PPzf  =  (A.22) 

This  can  be  viewed  as  minimizing  the  angle  between  the  waveform,  T,  and  the  row  space  of  the 
observations,  Z. 
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APPENDIX  B 

A  MAXIMIZATION  PROBLEM 


Consider  maximizing  the  determinant  ratio 


\X^AX\ 
T  \X»BX\ 


(B.l) 


over  all  full  rank  N  x  J  matrices  X  (J  <  rank(A)).  A  and  B  are  both  hermitian  positive  semi- 
definite  and  in  adoition  B  >  0.  Then 


\X^AX\  \X^B-'^f^AB-'^f^X\ 

\X»X\ 


"  |(X(X"X)-1/2)H(X(X«X)-1/2)| 

=  n  =  Yl  A,(5-M) 

*=1  fc=i 


where  S  runs  over  N  x  J  arrays  with  orthonormal  columns  (that  is,  satisfying  S^£  =  Ij).  The 
second  to  last  equality  is  well-known  (see,  for  example.  Bellman  [4],  p.  129).  This  inequality  can 
also  be  viewed  as  a  consequence  of  the  stronger  result  (Poincare  separation  theorem) 


found,  for  example,  in  Horn  and  Johnson  [5].  Letting  ej(-)  denote  the  N  x  J  matrix,  the  columns 
of  which  are  the  top  J  eigenvectors  of  its  argument,  the  maximum  occurs  when  X  —  ej(B“^A)r, 
where  T  is  an  arbitrary  nonsingular  J  x  J  array. 
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APPENDIX  C 
ARRAY  BEAMSHAPES 

The  angle  h  of  Equation  (7)  is  fundamental  in  the  geometrical  interpretation  of  signal  pro¬ 
cessing.  Mathematically,  6  is  proportional  to  the  geodesic  distance  between  the  two  array  responses 
and  ,  interpreted  as  points  in  projective  space.  Note  that  6  has  been  normalized  so 

that  6=1  corresponds  to  orthogonal  array  responses  in  the  noise-whitened  coordinates.  To  make 
physical  sense  of  this  definition,  consider  the  important  special  case  in  which  oc 

Some  new  notation  is  required.  For  the  most  part,  the  notation  in  this  section  is  not  used 
in  the  rest  of  the  report  and  any  conflicts  should  be  ignored.  Let  V{6,4>)  be  the  array  response 
parameterized  by  zenith  angle  (complement  of  elevation)  9  and  azimuth  angle  4>.  A  convenient 
coordinate  system  is  provided  by  the  unit  direction-of-arrival  vectors  u{9,  (^): 

(sin(5)cos(<^),sin(0)sin(^),cos(0))^ 

For  arrays  with  perfectly  matched  element  patterns  (a  common  modeling  assumption),  the  response 
of  the  fcth  element  can  be  written  {V  =  (ui,. .  .,vs)^) 

where  dk  is  the  three-dimensional  location  of  the  fcth  antenna  element  and  g{9,4>)  is  the  antenna 
pattern  common  to  all  elements.  Let 


UA  =  «(^l)^l)  - '“(^2,<^2)- 


Armed  with  this  notation,  one  has  the  approximation  (6  <C  1) 


where 


,7r. 


|V^"(^2,<^2)^^(^l,<^l)P 


1-  l|y(01,<^,)l|21|F(02>2)l|2 


a  1  -  -  dd^)'ttA 


N 


N-'^'£dkdl 

fc=i 
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and 


N 


k=l 


Thus  an  approximate  expression  for  b  becomes 
w  -  ddF)uA- 

When  UA  is  small,  this  expression  indicates  that  the  array  responses  at  angle  b  from  a  fixed  array 
response  V{62i<f>2)  fonn  an  ellipse  (planar  arrays)  in  (ux,Uy)  space.  Furthermore,  this  ellipse 
has  the  same  shape  no  matter  what  point  of  reference  V(02)<^2)  is  used.  By  extrapolating  this 
expression  out  to  6  =  1  [where  the  array  responses  and  V(92,<h)  would  be  orthogonal  if 

the  approximation  held]  one  has  the  definition  of  beamshape; 

l  =  ^ul(D-dd^)u^.  (C.l) 

This  equation  is  motivated  by  a  peak-to-nuU  definition  of  beamwidth.  As  mentioned  above,  the 
normalization  of  6  is  chosen  so  that  6=1  corresponds  to  a  null  (orthogonality)  in  the  pattern. 

For  planar  arrays,  Equation  (C.l)  describes  an  ellipse  in  {ux,Uy)  coordinates  with  major  and 
minor  semiaxes  given  by 

J  ■  -  dd^)  and  j  -  di’)- 

Equation  (C.l)  can  be  applied  to  arbitrary  array  geometries  to  yield  a  physical  beam.  For 
example,  a  completely  filled  (sampled  arbitrarily  finely)  line  array  of  length  A  has  a  beam  of  size 
y/Zfi  •  A/A.  A  filled  (including  interior)  circular  array  of  diameter  A  has  a  circular  beamshape 
of  radius  A/A.  In  all  cases,  the  beam  size  is  related  to  A/A  where  A  is  a  characteristic  size  of 
the  array.  Similar  beamshapes  can  be  motivated  by  Cramer-Rao  bounds  on  direction-finding  (see 
Delong  [6]). 
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APPENDIX  D 
PERFORMANCE 


It  is  difficult  to  characterize  analyticafiy  the  statistical  performance  of  the  tests  described 
above,  at  least  in  the  general  case.  However,  some  simple  coordinate  transformations  aUow  these 
statistics  to  be  cast  in  simpler  forms  that  depend  on  a  reduced  set  of  parameters.  These  normal 
forms  can  be  the  first  step  in  an  analytical  characterization  of  performance,  but  they  also  can  be 
used  as  the  basis  of  Monte-Carlo  evaluations.  For  the  most  part,  this  is  the  approach  taken  here. 

D.l  Normal  Forms 

Let  the  data  Z  have  covariance  R  ®  Ii,  as  above,  but  let  the  mean  of  Z  be  VAT  instead  of 
the  expected  VAT,  Here,  A  is  J  x  M  (instead  of  J  x  M)  and  T  is  M  x  L  (instead  of  M  x  L). 
This  models  a  mismatch  between  the  expected  signal  structure  and  the  true  structure.  V  has  the 
same  dimensionality  N  x  J,  but  represents  a  specific  value  of  V  below.  Choose  a  unitary  U  so  that 

where  ej,  is  the  column  vector  of  length  N  with  all  zero  entries  with  the  exception  of  a  single  one 
in  the  A:th  entry.  Define  2  UR~^I^Z.  Consider  the  partitioning 

(2p  2,)  =  2  ^  ^  j  =  {UR-^f^Zj,  vr-^'H^). 

The  covariances  of  these  new  variables  are 

cov{2)  =  In®  II 
cov(2p)  =  In®  Im 
cov(2,)  =  In®Il-m- 

It  is  convenient  to  introduce  some  notation  to  describe  their  means.  Define 
B  {V” R-'^Vfl'^AiTT^fl'^ 

C  B{TT^)~^I'^TP’^ 

D  B{TT^)-^I^TQ^ . 
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Then  the  means  of  the  new  variables  are 

E[2p]  =  EjC 
E[Z,]  =  EjD. 

D.2  Detection 

The  detection  statistic  can  be  expressed  equivalently  through  the  top  J  eigenvalues  of 

(D.l) 

This  form  is  convenient  for  characterizing  detection  performance.  Define  S  ^  .  In  terms  of 

the  new  random  variables,  one  can  write 

2^S-^2j,  =  =  Z^R-^Zj,. 

Thus  the  detection  statistic  can  be  expressed  in  a  normal  form  that  does  not  involve  the  scenario 
parameters  other  than  through  the  new  parameters  C  and  D. 

One  can  replace  2  with 


where  Um  and  Ul-m  unitaries.  This  causes  2p  to  be  replaced  by  2pUM  and  2q  by  ZqUi^M- 
The  detection  statistic  (depending  only  on  eigenvalues)  is  unchanged.  However,  the  p-mean  becomes 
EjCUm  while  the  q-mean  becomes  EjDUl-m-  Thus  detection  depends  only  on  the  values  of  CC^ 
and  DD^ .  These  values  are  determined  next. 

Define 


so  that 

CC^  = 

DD”  =  B{Im-TT^)B^  =bb" -CC^. 
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Note  that 


CC"  =  {V^R-^Vf<'^[AT{T^{TT^)-^T)7^A^]{V^R-^Vfl‘^  (D.2) 

<  {V^R-^VY/^{ATT^A^){V^R-^Vf^^  =  BB^  {D.3) 

with  equality  when  M  =  M.  and  T  =  T.  In  addition 

DD^  =  (V^R-^Vy^^[AT(lL-T^(TT^)-^T)T^A^](V^R-'^V)^^^ 

<  {V^ {ATT^ A^)iV^ =  BB^ . 

Detection  performance  depends  only  on  the  dimensional  parameters  (N,  L,  M,  J,  M),  the  nor¬ 
malized  signal  level  B,  and  the  mismatch  expressed,  for  example,  by  FF^. 

When  there  is  no  mismatch  (i.e.,  M  =  M  and  T  =  T),  one  can  show  that  the  detection 
statistic  depends  only  on  the  eigenvalues  of  CC^  (see  also  Kelly  and  Forsythe  [2]).  To  show  this, 
let  the  N  X  N  unitary  U  be  defined  so  that 


U  = 


Uj  0  \ 

0  In-j  )  ' 


Consider  changing  random  variables:  2  i->  1^2.  The  detection  statistic  is  based  on  the  top  J 
eigenvalues  of 


S-'^2p2^  ^US-^2p2flf^ 

and  hence  is  unchanged.  The  covariance  of  2  is  also  unchanged.  However,  the  p-mean  and  q-mean 
of  2  are  altered: 

E[U2p]  =  UEjC  =  EjUjC 
E[U2g]  =  UEjD  =  EjUjD  =  0. 

Thus  the  detection  statistic  depends  only  on  UjCC^U^  for  arbitrary  Uj,  and  hence  it  depends  on 
the  eigenvalues  of  CC^ ,  which  are,  from  Equation  (D.3),  the  eigenvalues  of  BB^ ,  namely,  those  of 

{V»R-^V){Arr^A^).  (D.4) 

WTien  J  =  1,  this  is  a  scalar  called  the  ideal  total  array  signal-to-interference  ratio  (TASIR). 
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D.3  Special  Case 

For  an  application  such  as  frequency  feature  processing  (FFP),  as  discussed  in  the  body  of 
the  text,  consider  the  special  case  when  J  =  I,  M  —  card(C),  M  ~  L,  T  is  unitary,  and  T  =  7c, 
the  array  formed  from  a  subset  C  of  rows  from  T.  Note  that  card(-)  counts  the  number  of  elements 
in  a  set.  The  rows  of  T  can  represent  any  orthouormal  signal  basis  (e.g.,  Fourier  components  or 
wavelets)  with  the  signal  of  interest  (SOI)  concentrated  in  a  particular  subspace  (i.e.,  subband  or 
time- frequency  cell). 

Let  Ap  denote  the  row  vector  of  amplitudes  corresponding  to  the  rows  indexed  by  C:  Ap  = 

, . . .  ,Ajj^).  Similarly  let  Aq  be  the  row  vector,  the  entries  of  which  come  from  the  complement 
of  C;  Aq  —  5  •  •  • , M )■  One  has 

E[2p]  =  eiC 

Y\Zq]  =  ejD 

with 

C  =  {V^R-^Vfl^Ap 
D  =  {V^R-^V)^^^Aq. 

Then 

CC"  =  (V"i2-'‘V)J3l>ljP  =  (v"7e-'V)pp|p  =  TASIRp 

i€C 

DD^  =  =  (v"7?~^V)|M,iP  =  TASIR, 

iic 

where  TASIRp  (TASIR,)  is  the  ideal  total  array  signaJ-to-interference  ratio  for  the  p-data  (q-data). 
From  Equation  (D.4),  TASIR  =  TASIRp  -I-  TASIR,. 

The  signal  powers  in  the  p-data  and  q-data  can  be  identified  with  the  averages  ll>tp|p/L  and 
ll^lgll^/Zr.  If  V  has  unit  norm  and  the  noise  floor  has  identity  covariance,  these  averages  are  called, 
more  precisely,  array  signal-to-noise  ratios  (ASNR)  (see  Section  4.1  for  definitions;  note  that  the 
noise  floor  covariance  is  denoted  there). 

D.3.1  Detection 

The  detection  statistics  shown  in  the  figures  are  essentially  constant  when  the  ideal  TASIR 
of  the  SOI  is  large.  A  heuristic  argument  can  be  given  to  approximate  this  constant  value. 
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Write  the  detection  statistic  as 

The  means  of  Zp  and  2,  can  be  written 

E[Zp]  =  e^iV^R-^V)^/\Ai,,...,Ai^)  =  er(V^R-^V)^/^Ap 
E[Z,]  =  eiiV^R-^V)^^\Ak,,...,Ak,_J  =  e,{V«R-^VYfU,. 

When 

a  M-\V^ R~^V)\\Apf  >  1, 

one  has 

M~^ZpZ^  ^aeie^. 

Unfortunately,  a  similar  approximation  is  not  always  appropriate  for  S.  Define 
f3  (L  -  M)-\V^R-^V)\\Agf. 

Let  Sp  be  the  random  matrix  obtained  by  replacing  {V^ R~^Vy^^Aki  P<>  drawn  from  a  set  of 
independent  identically  distributed  complex  Gaussian  random  scalars  with  zero  means  and  variances 
/9.  Then 


(L  -  M)-‘(V''fl-'V)||Af  =  (i  -  M)-'EE  IwPl  »  (i  -  M)-'  Z,  IwP 

I  I 

because,  when  L  —  M  >  1,  the  right-hand  side  has  a  complex  distribution  that  is  tightly 
distributed  about  its  mean.  Then,  because  the  distribution  of  S  depends  only  on  V)||>l5lP, 

not  on  the  specific  vector  Aq, 

{L  -  M)-^S  «  (I  -  M)-^Sp  w  (L  -  M)-^E[5p]  =  /jv  +  /Sejef . 


39 


Finally, 


«  qA^[(/a,  +  ^eief  r'ejef ]  =  a/(l  +  H). 

In  the  important  case  when  Af  =  1,  one  can  write,  when  X  »  1, 

TASIRp  ^  TASIRp  ,TASIRp 

a/(l  +  /^)  -  1  ^  TASIR,/(X  -  1)  ~  -^X  +  TASIR,  ^TASIR, 

where  TASIRp  and  TASIR,  are  the  total  ASIRs  in  the  p-data  (which  is  now  a  single  sample)  and 
q-data.  The  final  approximation  holds  when  TASIR,  >  X. 

D.3.2  False  Alarm  Statistic 

This  section  discusses  the  false  alarm  statistic  in  the  special  case  J  =  1  =  M.  Aside  from 
the  first  few  paragraphs,  it  is  self-contained;  the  notation  used  here  is  not  necessarily  related  to  the 
notation  in  the  rest  of  this  report. 

The  false  alarm  statistic  can  be  expressed  as(J=l  =  M) 

prob{X2^5-'2p  >r}. 

The  reason  for  introducing  the  factor  X  will  be  apparent  below.  Recall  that  Zp  and  2q  (and 
hence  S  =  ZqZ^)  are  independent.  Condition  on  the  value  of  Zp  and  choose  a  unitary  U  so  that 
UZp<x  e\.  Then 

ZfS-^Zp  =  {UZpfiUSU^)-\UZp)  =  \\Zpf{USU^)^^ 

where,  in  general,  denotes  the  (t,j)th  element  of  the  inverse  of  X.  But  USU^  and  S  have  the 
same  distribution  (denoted  USU^  ~  «S)  because 

cov{UZq)  =  cov(2,)  =  X;v  ®  Xl-i- 


Thus 


{USU^f^  ~  ~ 

X\L  -  N) 


(D.5) 
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where  x^iL  —  N)  is  a  complex  random  variable  of  degree  L—N.  Recall  that,  in  general,  x^(^)  has 
the  distribution  of  the  sum  H'^fclP  where  the  {zk}  are  independent  complex  circular  Gaussians 

of  zero  mean  and  unit  complex  variance.  The  last  equivalence  in  Equation  (D.5)  is  a  specific  case  of 
a  well-known  property  of  Wishart  distributions  (see,  for  example.  Appendix  1  of  Kelly  and  Forsythe 
[2]).  Summarizing  thus  far,  conditioned  on  one  has 


Z^S-^Zp  ~ 


xHL-NY 


Because  the  denominator  does  not  depend  on  the  conditioning,  one  can  write 


Z^S~^Zp  ~ 


X^L-NY 


Recall  that  the  incomplete  beta  function  is  defined 


Jo 


so  that  the  complete  beta  comes 


B{n,m)  jB(n,m,  1) 


r(n)r(m) 

r(n-t-m)‘ 


Integrating  by  parts,  one  has 

B^n,n,,a)  =  1.™-  I,,.). 

n  n 

Iterating  this  leads  to 


B(n,m,a)  =  X  n-^n  +  k)'  “  ‘ 


SO  that 


(D.6) 

(D.7) 
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The  left-hand  side  of  the  above  equation  is  the  cumulative  distribution  function  (CDF)  of  a  complex 
I3{n,m)  random  variable  (see  Kelly  and  Forsythe  [2]  for  a  thorough  discussion  of  the  complex 
random  variables  introduced  here).  This  random  variable  can  be  expressed  in  terms  of  a  complex 
random  variable  F{Tn,n)  which,  in  turn,  is  the  ratio  of  two  independent,  mean  zero,  complex 
random  variables: 


/3(n,  m)  ~ 


1 

1  -t-  F{m,n) 


F{m,n)  ~ 


X\rn] 

X^in)  ■ 


Consider  the  false  alarm  probability  for  the  detection  statistic  ( J  =  1  =  M) 
prob{Tf’(iV,X  -  N)>T}  =  prob{^(X  -N,N)< 


Using  Equation  (D.7)  one  has,  for  large  L  (which  is  the  case  of  interest  in  the  applications)  the 
limiting  behavior 


N—i  rpic 

Hm  prob{X  F(iV,  X  -  iV)  >  T}  =  e"^  V 


(D.8) 


This  expression  provides  the  large  X  false  alarm  behavior  shown  in  the  figures. 

The  false  alarm  probability  can  also  be  evaluated  by  employing  a  ChernofF  bound.  To  do  this, 
consider  the  moment  generating  function  for  the  log  of  the  beta  random  variable: 

EfgAln/?(n,m)j  ^  B{n,m)-^  r  <^+"-^(1  -  dt  =  + 

Jo 

_  r(A  +  n)  r(n  -f  m)  _  n  +  k 
r(A-fn-l-m)  r(n)  j"o  k  +  X’ 

Let  Q  In  Then 

prob{/?(n,m)  <  =  prob{ln /3(n,m)  <  Q}. 
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Define 


TTl  — 1  ,  » 

/(A)  H  "t  "V 

L_A  +  K  -  A 


Chernoflf’s  bound  on  the  distribution  tail  is  given  by 


prob{ln;3  <Q)<  inf  =  inf  /(A). 


The  derivative  of  /(A)  is  ^ven  by 


m— 1 

/(^)  =  [^?  +  E 


1 


fc=0 


n  +  fc  —  A 


]/(A). 


It  is  clear  that  the  stationary  points  are  relative  minima.  One  can  thus  express  the  Chernoff  bound 
parametrically  as 


Q{X) 

m 

T 


m— 1 


-E 


1 

n  +  k  ~  X 


m— I 

eAQ(A) 

fc=0 


n  +  k 
n  +  -  A 


-  1. 


(D.9) 

(D.IO) 

(D.ll) 


A  limiting  form  of  the  Chernoff  bound  is  relevant  for  this  report.  Let  m  =  N  and  n  =  L  -  N  for 
the  application  to  FFP.  Further,  replace  T  with  TjL  where  £  >■  1.  Then 


prob{/?(£ -iV,JV)< 


1 


1  +  T/L 


N-\ 

}  <  n 

k=0 


L-N  +  k 
L-N  +  k-X 


and 


iV-l 


L-N  +  k 
L- A  +  fc-A 


e 


-r 
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when  T  >  N,  because  Amin  ~  - and  Q(A)  w  -T/L  (really,  AminQ(Amm)  -*  N  -  T  and 

Amin/ ^  1- 

Figure  D-1  compares  the  asymptotic  expressions  and  bounds  discussed  above  to  the  true  false 
alarm  probability  in  one  example.  For  this  example,  L  =  1000  and  iV  =  4.  The  Chernoff  bound 
is  given  by  the  parametric  expression  above.  The  asymptotics  of  this  expression  are  also  shown. 
It  is  clear  that  the  asymptotics  of  the  incomplete  beta  [Equation  (D.8)]  are  quite  accurate;  these 
asymptotics  are  used  for  Figure  3. 

D.4  Copy- 

Copy  can  be  understood  in  terms  of  signal-to-interference  ratios  (SIR).  This  section  considers 
very  general  notions  of  SIRs  based  on  the  estimate  of  the  parameters  A.  These  SIRs  are  specialized 
to  FFP  applications  in  a  subsection. 


230287-1 S 


Figure  D-1.  Comparison  of  Incomplete  Beta  function,  bounds,  and  asymptotics. 


Some  additional  notation  is  required.  Define  the  JK  column  vector  Vec(A')  of  the  J  X  K 
matrix  A'  so  that 

Xjk  =  (Vec(X))(j_i)/c+i.. 

Some  properties  of  Vec  are  easily  verified  for  arbitrary  matrices  (with  conformal  dimensions) 
X,  y,  c,  D-. 

1.  Vec(A')^Vec(y)  =  tr(A'"Y) 

2.  {C  ®  D)Vec(X)  =  YeciCXD'^) 

3.  Vec(y)«(C  ®  D)Yec(X)  =  tr(y"CXD^)  . 

Recall  that  Kronecker  products  are  defined  in  Section  A.l.  The  ML  estimate  A  of  the  signal 
amplitudes  A  is  ^ven  by 

A  =  W^ZT^iTT^)-^  =  {W®  T^(fT^P‘)"Vec(Z).  (D.12) 

Thus  one  can  view  both  W  and  T  as  weights.  To  form  W,  one  needs  an  estimate  V  of  V,  provided 

by 

V  =  ejiZpZ^S-^)e 

where  ejr(-)  denotes  the  N  x  J  array  formed  by  the  top  J  eigenvectors  of  its  argument.  Only  the 
column  span  of  ej  is  unique.  The  ambiguity  is  expressed  by  the  post-factor  0,  which  changes  bzises 
for  the  column  space.  Because,  in  general, 

w  = 

one  has 

W  =  ej(S-'ZpZf  )0. 

Again,  the  postfactor  0  (different  each  time)  expresses  nonuniqueness.  Define 
W  =  ej{S-^ZpZ^)Q. 
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Ideally,  one  would  characterize  the  estimate  A  by  providing  distributional  information,  or 
failing  that,  first  and  second  moment  information.  A  less  ambitious,  but  physically  reasonable 
approach  is  taken  here.  In  this  approach  W  (and  hence  W)  is  assumed  constant  in  the  calculation 
of  a  SIR  based  on  first  and  second  moments  of  A.  Then  the  statistical  behavior  of  W  is  introduced 
to  characterize  this  SIR  as  a  random  variable.  This  makes  sense  if,  for  example,  the  data  used 
to  determine  W  (and  T)  and  the  data  used  in  Equation  (D.12)  are  independent.  Then,  using  the 
properties  of  Vec  described  above  and  the  properties  of  covariance  described  in  Section  2.1,  one  can 
write 

covv^(A)  =  (W^RW)  ®  iTT»)~^  =  (>V">V)  0  (Tr")"\ 

In  addition, 

Vec(Evv[i])  =  Vec[W^VA{TT^){TT^)-^]  =  Vec[yV^  EjBT{TT^)-'^^^] 

=  (W  ®  (TT")'^^^)^  Vec[EjBT], 

These  “mean”  and  “covariance”  statistics  are  based  on  a  fixed  value  of  IT;  hence  the  W  subscripts. 
One  can  define  a  random  variable  called  the  achieved  TASIR  as  follows: 

"  Vec(E,vlil)''covw(Ar‘Vec(Ewl4)). 

Upon  using  the  properties  of  the  Kronecker  product  and  of  Vec,  this  becomes 

VeciEjBTf  (Vy(>V^)V)-^>V^  ®  /a/)  Vec(£j6r) 

=  Yec{Bf  [{E^WiW^Wy'yV^Ej)  ®  {TT^f]  Vec(B) 

=  tr[5^i;^>v()v^>v)-‘>v^EjBrr^) 

=  tr[i;5^>V(>V"H’)-^)V^£;jCC"] 

where  C  is  defined  in  Section  D.l. 

Note  that  depends  only  on  the  column  space  of  W,  because  it  is  the  or¬ 

thogonal  projector  onto  that  subspace.  Furthermore,  the  argument  of  Section  D.2  shows  that  the 
distribution  of  this  projector  only  depends  on  the  dimensional  parameters  N,  M,  J,  L,  and  on  the 
signal  levels  expressed  through  CC^  and  DD^ .  Because  the  projector  is  the  only  random  variable 
in  the  achieved  TASIR  expression  above,  this  is  the  desired  normal  form  for  copy  performance. 
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More  generally,  one  is  interested  in  the  SIR  based  on 


where  W  (and  T)  are  determined  as  above,  but  the  M  x  L  array  T  “selects”  a  different  part  of  A 
(see  special  Ccise  below).  Defining  C  so  that 

cc^  =  (v"i^-^v)^/2[.Ar(f"(ff")-^f’)r^>l"](v"R-^v)^/^  (d.i3) 

the  calculation  above  yields  an  achieved  TASIR  given  by 

tT[E^WiW^W)-^W^  EjCC^]. 

As  before,  W  depends  only  on  the  dimensional  parameters,  CC^,  and  DD^ . 

The  reason  for  the  more  general  SIR  is  that  the  weights  W  and  T  may  be  determined  by  only 
part  of  the  signal  (the  spectral  lines  in  FFP  applications)  but  copy  performance  may  be  required 
on  a  different  part. 

D.4.1  Special  Case 

Recalling  the  expression  for  CC^  above,  when  J  =  1,  the  achieved  TASIR  becomes 

(ef  >V(>V^>V)-^W"ei)  (V"R-‘V)  (>irf^(f’f’^)-’f  r">l") 

=  (efW(>V">V)-^>V"ei)(V^i2-^V) 

k^C 

where  T  =  (see  Section  D.3  for  notation).  The  first  factor  expresses  the  loss  due  to  sampling 
(including  suppression  effects  due  to  mismatch,  as  expressed  by  DD^).  The  second  factor  (if  V  has 
unity  norm)  expresses  the  loss  of  SIR  due  to  the  spatial  distribution  of  interference,  and  the  last 
factor  is  the  signed  energy  in  the  subband  expressed  by  the  set  C.  It  is  important  to  note  that  the 
loss  factor  is  common  to  all  SIR  expressions,  no  matter  what  part  C  of  the  signal  is  used  to  assess 
copy.  In  the  particular  case  when  C  =  C,  the  above  can  be  written 

(ef  W(yv"VV)-^>V"ei)  TASIRp. 

Of  course,  the  achieved  ASIR,  which  is  defined  on  a  per  look  basis,  is  always  smaller  by  a  factor  of 

L. 
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